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Abstract 

The presently experimentally interesting form factor for the anomalous pro- 
cess 7 — > 7r + 7r°7r - is calculated as the quark "box" -amplitude where the inter- 
mediate fermion loop is the one of constituent quarks with the pseudoscalar 
coupling to pions. This also corresponds to the form factor, in the lowest 
order in pion interactions, of the cr-model and of the chiral quark model. We 
give the analytic expression for the form factor in terms of an expansion in 
the pion momenta up to the order 0(p s ) relative to the soft point result, and 
also perform its exact numerical evaluation. We compare our predictions with 
those of the vector meson dominance and chiral perturbation theory, as well 
as with the scarce data available so far. 

PACS: 13.40.Gp; 14.40.Aq; 12.38.Lg; 24.85.+p 
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Typeset using REVTeX 



The analysis of the Abelian axial anomaly shows that the 7r° — > 77 amplitude is 
exactly T 2l (m n = 0) = e 2 iV c /(127r 2 / 7r ) in the chiral and soft limit of pions of vanishing 
mass m n . (Here, e is the proton charge, and N c = 3 is the number of quark colors.) 
This amplitude T 27 (m 7r = 0) is successfully reproduced also by the simple "free" quark 
loop (QL) calculation of the pseudoscalar- vector- vector (PVV) "triangle", provided one 
uses the quark-level Goldberger-Treiman (GT) relation g/M = l/f n to express the effective 
constituent quark mass M and quark-pion coupling strength g in terms of the pion decay 
constant /„. = 92.4 MeV. This calculation (essentially a la Steinberger ||) is the same as 
the lowest (one-loop) order calculation in the a-model which was constructed to realize 
current algebra explicitly. By "free" quarks we mean that there are no interactions between 
the effective constituent quarks in the loop, while they do couple to external fields, presently 
the photons and the pion 7r a . Our effective QL model Lagrangian is thus 

Ctf = V {ip- eQ4- M)V - ig VWaTa* + ... , (1) 

where Q = diag(Q u , Qd) = diag(|, — |) is the quark charge matrix, and r a are the Pauli 
SU (2)-isospin matrices acting on the quark iso-doublets ^ = (u, d) T . The ellipsis in C e g 
serve to remind us that Eq. (|I|) also represents the lowest order terms from the u-model 
Lagrangian which are pertinent for calculating photon-pion processes. The same holds 
for all chiral quark models (xQM) - considered in, e.g., Ref. - which contain quark- 
meson coupling M^(UPl + WPr)^ with Pl,r = (1± 7s)/2. Namely, expanding U®> = 
exp[(— )iii a Ta/ fn\ to the lowest order in 7r a and invoking the GT relation, again returns the 
QL model Lagrangian ([[]). 

This simple QL model (and hence also the lowest order xQM and the a-model) provides 
an analytic expression {e.g., see Ref. 0) for the amplitude T 27 (m 7r ) also for m n > (but 
restricted to m n < 2M, which anyway must hold for the light, pseudo-Goldstone pion), 
namely 
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Adler et ai, Terentev, and Aviv and Zee [^] proved that the amplitude for the anomalous 
processes of the type 7 — ^ 37r is related to T 27 (0) and is given by 

^(0,0,0) = ^T 27 (0) = (3) 

The arguments of the anomalous amplitude (H), namely the momenta {^1,^2,^3} of the 
three pions {vr + , tt°, vr~}, are all set to zero, because Eq. (^) is also a soft limit and chiral 
limit result, giving the form factor F^ n (pi,p 2 ,P3) at the soft point. 

In the QL model, the amplitude (|3[) is obtained by calculating the "box" graph, Fig. 1. 
This is not surprising, as the anomalous "box" amplitude (0) was already obtained analyti- 
cally and exactly by Alkofer and Roberts [0 in the so-called Schwinger-Dyson (SD) Ansatz 
approach, which is more general. In this approach, one employs an Ansatz for a dynamically 
dressed quark propagator, characterized by a momentum-dependent quark mass function 
M(k 2 ), and by the related momentum-dependent quark- ant iquark pseudoscalar pion bound 
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state Bethe-Salpeter vertices IV instead of our g^T a quark-pion Yukawa couplings. The 
present free QL model, with its constant quark-pion coupling strength g and the free quark 
propagators S(k) = 2(7 • k + M)/{k 2 — M 2 ) containing constant effective constituent mass 
M, can therefore be considered as a special case of the SD Ansatz approach, i.e., the simplest 
possible such Ansatz. 

Let us stress that for the processes of the type 7 — > 3tc, going beyond the soft limit 
is much more important than for the process 7r° — ► 27 where the amplitude T 2 " 1 obtained 
in the chiral and soft limit is an excellent approximation for the realistic 7r° — > 77 decay. 
By contrast, the current TJNAF measurement of the 77r + — ► 7r + 7r° process ||, as well as 
already published || and still planned Primakoff measurements at CERN ||10|| , involve so 
large values of energy and momentum transfer, that departures from the soft-point result 
(0) may well be significant. 

Nevertheless, to the best of our knowledge, no extension of Eq. (^) beyond the chiral 
and soft limit has been given in the QL (%QM, a-model) approach so far, neither through a 
numerical calculation, nor as an analytic result that would be the analogy of the amplitude 
(0) for the 77-decay of the massive pion. (In the SD approach J7[, the box graph beyond 
the chiral and soft limit was addressed, but only numerically.) To fill this gap that seems 
to exist in the literature even at the lowest non-trivial order, we provide such an analytic 
expression in the form of an expansion in powers of the pion momenta up to the eighth 
order, in addition to the numerical evaluation of the 7 — > 3n form factor accurate to all 
orders. 

To compute the 7 — > 3it amplitude we use Eq. ([!]) in terms of the physical fields ir^ = 
(tti =F27T2)/\/2 and 7r° = 7r 3 , so that 71777 = ^/2(n + r + + tt~t_) +7r°r 3 where r± = (17 ±zr 2 )/2. 
There are six different contributing graphs, obtained from Fig. 1 by the permutations of the 
vertices of the three different pions. The momenta flowing through the four sections of the 
quark loop are conveniently given by various combinations of the symbols a, /3, 7 = +, 0, — 
in k af3l = k + api + /3p 2 + 7P3. Then, the anomalous vertex Vp coupling 7r + 7r°7r~ to 7, viz. 
the scalar form factor F^ 7T (p 1 ,p 2 ,p 3 ) associated to it, is in the present approach calculated 
through the six VPPP box graphs as 

Vf, = ie^pv p\f % vl F^ n (pi,p l ,p 3 ) 

= itfiupa P1P2P3 fj*(PhP2,P3) + (permutations of n + , it , tt~) , (4) 
where the contribution of the first diagram (Fig. 1), ie^pu PiP^Psf^ : {pi,P2,Pz)i is given by 
f d^k 

- J Tr {eQ7^ S(k) iV2 g>y 5 T + S(k„ 00 ) ig<y 5 r 3 S(k__ ) iV2 g>y 5 T_ S(k — ) j . (5) 

The color and isospin traces are N c and Ttj(Qt + t 3 t^) = —2/3, respectively. The 
Dirac trace is AiMe^pa P1P2P3 , leading to the partial amplitude /~ 7r (pi,P2,P3) = 
(+2/3)(eN c g 3 M/2ir 2 )I(p u p 2 ,p 3 ), where 

I(Pi,P2,P3) = ~^2 J (p _ M 2 )(k 2 00 - M 2 )(k 2 __ - M 2 )(F__ - M 2 ) ' 

After combining the four propagator denominators by the Feynman trick, and shifting k^ 
by pi(l — x 3 ) +p 2 (l — x 2 ) +^3(1 — xx), the integral over the loop momentum k becomes 
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where 



R(,Pi,P2,P3) 2 = M 2 -plxi(l - xi) -p\x 2 (\ - x 2 ) -p\x 3 (l - x 3 ) 

-2p 2 ■ p 3 x 2 (l - xi) - 2pi ■ p 2 x 3 (l - x 2 ) - 2pi • p 3 x 3 {l - xi) . 



(8) 



Adding up the contributions from the remaining 5 diagrams yields the total amplitude 
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Since the integrals such as I(pi,P2,P3) are symmetric in the interchange of the first and the 
third argument, the two curly brackets in Eq. (|Sp are equal to each other. Therefore, the sum 
of the three w-quark "box" loop diagrams, as well as the sum of the three <i-quark ones, give 
contributions to the amplitude © which are separately symmetric under p\ <-> p 2 <-> p 3 . 
For that reason, a calculation that would employ the integer charge iso-doublet nucleons 
in the loops instead of the ones with quarks, would - up to the fermion mass values, M 
vs. M nuc i eon - give the same amplitude in spite of the three neutron graphs dropping out 
because of the vanishing neutron charge. Also, the interplay of the charge-isospin factors 
compensates the quark color factor N c = 3 as in the 7r° — > 77 amplitude T^ 7 , although the 
quark charges enter in F 3n linearly, and not quadratically. 

The integrals appearing in Eq. (|9[), exemplified above by the one defined by Eqs. (§D~(§|), 
are calculated in two ways. Firstly, we calculate them numerically, by the Gauss-Kronrod 
method. Secondly, we obtain the analytic expressions for F^ [Eqs. fll0D,(|l2T), (|13D and (|14D 
below] by expanding the integrands [exemplified by the one in the integral @-@] in the 
series of the scalar products of the pion momenta, Pi ■ pj = 1,2,3). In this case, the 
analytic integration over Xi,x 2 and x 3 finally yields [to order 0(p 4 )} 
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After using the GT relation, Eq. ([H]) returns at the soft point (pi = 0) the axial 
anomaly result (Q). We introduce the form factor normalized to the anomaly amplitude @, 
F^ n (pi,P2,P3) = F^ w {pi,p2,P3)/F^ K (0, 0, 0). It is also convenient to re-express the scalar 
products Pi ■ pj through the Mandelstam variables as defined by Ref. ||: s = (pi + P2) 2 , 
t' = (P2 + P3) 2 , u = (pi + P3) 2 , while t = p\ serves as the measure of virtuality of the third 
pion which may be off shell. The photon momentum is q = p\ + P2 + P3- 

We obtained the expansion for F^ n {p\,p2,pz) to the order 0(p 8 ) [relative to the anomaly 
result @], but give Eq. ( |I~0"|) only to the order 0(p 4 ) for brevity as the C(p 8 )-expansion for 
general p^ is very lengthy. It will be given elsewhere. However, below we do give F^ n to the 
order 0(p 8 ) for the simpler special cases which are important for comparing our predictions 
with the experiments [pjpjJlOl. Namely, one can take the photon to be on shell in all three 
pertinent experiments PJ5,lQ"f, in which at least two pion momenta, those of ir + and 7r°, are 
also on the mass-shell. We thus set q 2 = and p\ = p\ = m 2 ,, whereby 

s + t' + u = p 2 1 +p 2 2 +p 2 3 + q 2 = 2m 2 n + t. (11) 

In the Primakoff measurements [PJlCfl, including the one [[| providing the only existing 
data point so far, the third pion is also on shell. Hence, t = m 2 , in which case we predict 

P3*( s j/) = (i+ m * + m * + imm l + 193 ™* 
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+ 1 -— (s 4 + 2sH' + 3s 2 t' 2 + 2s t' 3 + t' 4 ) + 0(p 10 ) . (12) 



At CEBAF ||, one takes data near t ~ —ml, for which our 0(p ) expansion yields 
p*( a ^ = (i + m * + m ^ + 13 m ^ + m ^ ^ _ m ^ ^ 8 '"- 
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where the s t' symmetry is lost due to the virtuality t = —p\ = —p\ = —ml of the third 
pion. 
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Our momentum expansions flT2|) and (|T3|) show clearly that the main contribution to 
the term linear in s and t' (dominating the s, t'-dependence close to the soft limit), comes 
from 0{p A ) and not 0(p 2 ). This happens since the C(p 2 )-terms, due to the constraint (PD, 
contribute only to the part independent of s and t'. 

Note that in the on-shell case (0), the finite pion mass m n causes a larger upward shift 
than in the off-shell case ([13]). For chiral pions (m n = 0) and real photons the condition 
(O) becomes s + t! + u = t. For this case, but for general t, the amplitude (0) becomes 



F' (s,t) l+o« + oni#d + i/ni/fi + con n Jrt ^ + 7 J« + 7 Jrt ( S + * ) 




945 M 8 

1890 M« V + 2s¥ + 3s¥2 + 2s f ' 3 + *") + °^°) > ( 14 ) 
showing that the s <-> t' symmetry is restored in the chiral limit. The massless pion amplitude 



([14]) is smaller than the one with m n = 138.5 MeV by, typically, 4% when M = m p /2 = 385 
MeV, by some (depending on s) 6% when M = 330 MeV, by more than 10% when M = 250 
MeV, etc. It is interesting that for small s and t', the chiral-limit F^ n can fall slightly below 
its soft point value (|3|) when t < 0. 

In the CEBAF experiment ||, s will vary more than t'. For the t'-range relevant at 
CEBAF, the t'-dependence of F^ n anyway turns out to be rather weak. For example, 
suppose one plots (not done here to avoid overcrowding our figures) the t = — m 2 form 
factor (]T3| ) as a function of s for several values of t' varying from t' = —m 2 to t' = —8m 2 . 
(Take M = 330 MeV for definiteness.) One would thus get, across the whole s-range relevant 
at CEBAF, a narrow strip of tightly spaced curves, where the curve depicting F^(s) for 
t! = —8m 2 , would differ by just 2% to 3% from the curve with t' = — m 2 . In Figs. 2 and 
3, we therefore show the dependence of F^ n (for various cases) on the variable s, with t' 
fixed. On the other hand, since t changes sign when we put also the third pion off-shell as 
at CEBAF ||, the amplitude is more sensitive to this change. 

In Fig. 2 we give our results at t — — m 2 (the t- value most relevant at CEBAF ||) for the 
constituent masses M = 385 MeV, 330 MeV, 300 MeV and M = 250 MeV. In fact, we show 
a pair of curves for each of these masses: one of the curves is obtained by exact numerical 
evaluation of the box amplitude, and the other corresponds to our series expansion to the 
order 0(p 8 ). At the lowest depicted s, the form factor obtained by the accurate numerical 
calculation is (for each M) slightly below the corresponding series expansion approximating 
it, but exceeds it eventually as s grows. The convergence of the expansion is very satisfactory 
on the whole, since for s<llm 2 , the agreement between the two ways of calculation is very 
good for all these values of M. For s > 11m 2 , the exception are only the cases with 
unrealistically small M, such as M = 250 MeV. For M = 250 MeV, we plot the curves 
up to s — s tT w 13.03m 2 only. Namely, when s reaches s tr = (2M) 2 , i.e., the threshold 
for production of an on-shell quark- ant iquark pair, the QL approach becomes inadequate, 
because the amplitude starts being dominated by this threshold which is not physical but an 
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artifact of the model (|l|). [Concerning the accuracy of the computations close to a threshold: 
the difficulty of numerical integration starts increasing gradually as one gets closer than m 2 
to the threshold, while the accuracy of the 0(p 8 ) expansion starts failing before that.] 

Nevertheless, for the values of s < 16m^ accessible at CEBAF ||, such a threshold 
cannot be reached unless M < 2m 7r = 277 MeV. Such values are, however, too low to serve 
as the constituent quark masses, which cannot be much lighter than M « M nuc i eon /3 ~ 313 
MeV. 

In Fig. 3 we compare our (numerically obtained) on-shell predictions for various M with 
the only existing experimental point so far |J, and with the predictions of chiral perturbation 
theory (xPT) [|llj (in Holstein's ]12| renormalization convention - i.e., we take his [12] Eq. 
(10) for the xPT form factor), of vector meson dominance (VMD) fl3| [Holstein's [0] Eq. 
(9)], and of VMD with the final pion rescattering included fT2|j . 

We conclude that our results agree rather well with VMD for constituent quark masses 
300 MeV < M < 330 MeV when s > 8m 2 . Going down in s, already at s ~ 8m^ we 
agree rather well, but for somewhat higher M (330 MeV ^M<m p /2), with both VMD 
and xPT. For all s-values shown, our predictions get somewhat closer to those of xPT 
when the ratio m\jM 2 gets smaller, i.e., for the largest considered constituent quark mass, 
M = vripjl = 385 MeV. Since xPT results in the weakest momentum dependence, its 
predictions for the largest values of s accessible at CEBAF are significantly different from 
both VMD and our QL (xQM, cr-model) approach. The CEBAF experiment should thus be 
able to distinguish between these various physical mechanisms in this range of momenta. On 
the other hand, all these approaches agree reasonably well for the lowest of s-values accessible 
at CEBAF. In particular, our approach agrees with VMD and xPT that the existing data 
point [0 is probably an overestimate, as we can fit it well (e.g., with the M = 250 MeV 
curve in Fig. 3) only for unrealistically low values of M. 
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FIGURE CAPTIONS 



Fig. 1: One of the six box diagrams for the process 7 — > 7r + 7r 7r~. The position of the u 
and d quark flavors on the internal lines, as well as Q u or quark charges in the 
quark-photon vertex, varies from graph to graph, depending on the position of the 
quark-pion vertices. 

Fig. 2: Our numerically obtained 7 — * 3n form factors are compared with the corresponding 
O (p 8 ) -expansions (O) for various values of M and m^ = 138.5 MeV. The correspond- 



ing pair of curves for M = 385 MeV is denoted by the solid lines, by the short-dashed 
lines for M = 330 MeV, by the long-dashed lines for M = 300 MeV, and by the very 
short-dashed lines (only up to s tT = 13.03 m 2 , denoted by x's) for M = 250 MeV. The 
curves resulting from numerical integration start slightly below their respective power 
series at lowest values of the s-variable, but then exceed the latter for sufficiently high 
s. All curves pertain to the off-shell case t = — m 2 (dominant at CEBAF ||, but 
not in the Serpukhov experiment which provided the displayed data point, where 
t = ml). The remaining variable if is also set to if = —vn 2 ^. 

Fig. 3: The F^(s,t' = -ml) predicted by VMD (solid curve), VMD with final pion 



interactions 0] (dotted curve), and |]ll~| , |T2"f (dash-dotted curve), are compared 
with our F^(s, if = —m^) obtained by numerical integration for M = 330 MeV (short- 
dashed curve), M = 300 MeV (long-dashed curve), and M = 250 MeV (the topmost, 
very short-dashed curve). Same as the displayed data point |J, all curves pertain to 
all three pions on shell, so that t = m^ = (138.5 MeV) 2 . The remaining free variable, 
t\ is set to t' = —ml for all curves for definiteness. 
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